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ABSTRACT 

Moderate resolution, near-IR spectroscopy of MS1512-cB58 is presented, 
obtained during commissioning of the the Near IR Spectrometer (NIRSPEC) 
on the Keck II telescope. The strong lensing of this z = 2.72 galaxy by the 
foreground cluster MS1512+36 makes it the best candidate for detailed study of 
the rest-frame optical properties of Lyman Break Galaxies. In eighty minutes of 
on-source integration, we have detected Ha, [N ll]A6583,6548A, [O l]A6300A, 
He I A5876A, [O lll]A5007,4959A, E/3, H7, [O ll]A3727, and a strong continuum 
signal in the range 1.29-2.46 /xm. 

A redshift of z = 2.7290 ± 0.0007 is inferred from the emission lines, in 
contrast to the z = 2.7233 calculated from UV observations of interstellar 
absorption lines. Using the Balmer line ratios, we find an extinction of E(B- 
V)=0.27. Using the line strengths, we infer an SFR= 620 ± 18 M /yr(H o =75, 
qo=0.1, A = 0), a factor of 2 higher than that measured from narrow-band 
imaging observations of the galaxy, but a factor of almost 4 lower than the SFR 
inferred from the UV continuum luminosity. The width of the Balmer lines 
yields a mass of M V i r = 1.2 x 10 10 M . We find that the oxygen abundance is 
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1/3 solar, in good agreement with other estimates of the metallicity. However, 
we infer a high nitrogen abundance, which may argue for the presence of an 
older stellar population. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
individual: MS1512-cB58 



1. Introduction 

The largest sample of high redshift galaxies was selected from observations of the 
extinction of rest-frame far-ultraviolet light by intrinsic and intergalactic absorption. These 
Lyman Break galaxies (hereafter LBGs; Steidel et al. 1996) may, in principle, be the 
tracers of the global star formation history of the universe (Madau et al. 1998). However, 
the unquantified effects of dust extinction present an obstacle to such interpretation. For 
example, corrections to the star formation rate (SFR) of LBGs have been suggested to be 
as large as factors of 2-10 (Pettini et al. 1998, hereafter P98; Trager et al. 1997). 

In this paper we present the first observations of an LBG with NIRSPEC, the near 
infrared spectrometer on the Keck II telescope. We will use the rest-frame optical spectrum 
to obtain a model independent measure of dust extinction, an estimate of the virial mass, 
and the star formation rate. Spectra such as this one also provide the best method for 
determining the metal abundance in LBGs. We will examine the evolutionary status of 
cB58, and its implications as a "typical" LBG. 

We have chosen the brightest LBG as our first target (MS1512-cB58; see Yee et 
al. 1996). MS1512-cB58 is a gravitationally lensed starburst galaxy at z=2.72. It was 
discovered serendipitously, within the field of the z=0.37 galaxy cluster MS1512+36, in 
spectra obtained for the CNOC survey (Yee, Ellingson, & Carlberg 1996). MS1512-cB58 
is clearly extended with the morphology of a lensed arc (Williams & Lewis, 1997; Seitz et 
al., 1998). It has an ultraviolet spectrum representative of LBGs but due to a factor of 
~ 30 magnification (Seitz et al. 1998) it is orders of magnitude brighter than any other 
object of its kind (V=20.6, K'=17.8; Ellingson et al. 1996, hereafter E96). Like most LBGs, 
its UV spectrum (devoid of strong, high ionization emission lines) rules out any non-stellar 
(AGN) contribution to its flux. Its apparent (lensed) star- formation rate is 2417 M /yr (H 
= 75 km s~ 1 Mpc _1 , qo = 0.1; with extinction correction) from the 1500A continuum 
(Pettini et al. 2000, hereafter P2000). Throughout this paper, unless otherwise noted, a 
cosmology of (H = 75 km s _1 Mpc~ 1 , q = 0.1, A = 0) is assumed. 
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2. Observations 

Near IR spectra were obtained using the NIRSPEC instrument on the 10-m W. M. 
Keck II telescope (McLean et al. 1998, 2000a). Spectra were obtained in the low resolution, 
long-slit mode, using a slitwidth of 3 pixels (0.57") at the (InSb) detector. Observations were 
made in the H and K bands, as well as the wavelength range between J and H (1.3-1.6 /im; 
hereafter N4 for the custom NIRSPEC-4 filter) The seeing was 0.4", 0.4", 0.6" for N4, H, 
K respectively. MS1512-cB58 is highly elongated with a major axis ~ 2" long, but a minor 
axis only a few tenths of an arcsecond wide. The slit was aligned with the long axis, while 
the unresolved (seeing broadened) width filled the slit. H-band spectra were obtained in 300 
second exposures separated by ~ 10" nodding along the slit; N4-band and K-band spectra 
were taken in 600 second integrations. Total integration times were 1200, 1800, and 1800 
seconds in N4, H,and K respectively. 

The data reduction procedure followed the steps outlined in McLean et al. (2000b). 
Wavelength calibration was good to 0.99A. The output pixel scale was 4.16A/pixel in H 
and K, and 2.08A/pix in N4. Unresolved Argon arc-lamp lines had FWHM= 13. 7A, for a 
final resolution of, for example, R= 1300 at 1.8 /im. Nodded frames were registered, and 
the final ID spectrum was extracted using a Gaussian optimal extraction filter. Extracted 
spectra were divided by the atmospheric absorption spectrum. We obtained this spectrum 
from similarly- reduced observation of a mostly featureless star (HR5569, an A2V for the 
K-band; PPM130690, a G star for the H-band; and HR5630, an F8V star for the N4-band), 
divided by a Kurucz (1993) model atmosphere. Emission-line equivalent widths (EW) were 
measured using Gaussian fits to the line after subtraction of a local continuum. The flux of 
the line is then the EW multiplied by the continuum flux density, which is well known from 
broad-band photometry (E96). MS1512-cB58 is a homogenous (featureless) lensed image of 
the highest surface brightness part of the source galaxy (Seitz et al. 1998); as such, our 
measurements of the EW can scale directly to the photometry, since light lost outside the 
slit is just a pure fraction of the total light. No real features fall outside the slit, as might 
be expected for an unlensed galaxy. 



3. Results 

Figures 1-3 show the complete 1.28-1.53 /im, 1.55-1.97 /im and 2.05-2.46 /im spectra 
of MS1512-cB58. Just as LBG spectra in the rest-frame UV are remarkably similar to 
nearby starbursts, so too the optical spectrum of MS1512-cB58 is very similar to local 
irregulars or late type spirals. Table [l] lists the wavelengths and fluxes of the measured 
emission lines. 
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The mean redshift calculated from the emission lines is 2.7290 ± 0.0007, in contrast 
to the z abs = 2.7233 ± 0014 reported in Yee et al. (1996) from the rest-frame UV 
interstellar absorption lines. P98 find that velocity shifts on the order of several hundred km 
s _1 between interstellar absorption lines and H II emission lines are typical in LBGs. P2000 
confirm that this shift is seen in MS1512-cB58 as well, finding z em = 2.7296 ±0.0012 in good 
agreement with our result. Since the measurement of the systemic redshift of MS1512-cB58 
has only recently been available, several observations have been made in recent years 
assuming z em = z a b s (e.g. Bechtold et al. 1997; Frayer et al. 1997;.Nakanishi et al. 1997). 
In most cases the observations are still valid despite the shift. 

P2000 find evidence in the UV spectrum for intervening absorption systems at 
z=0.829 and z=1.339. The optical spectrum shows no emission features associated with 
those systems. However, the absorbing galaxy might be several arcseconds away from 
MS1512-cB58, and so might not have fallen in the slit. 

3.1. Dust Extinction and Virial Mass 

The extinction correction due to dust is crucial to understanding the star formation 
activity occurring in the object. E96 infer E(B-V)=0.3 from broad band photometry 
compared to synthetic galaxy spectra, and the LMC extinction curve of Fitzpatrick (1986). 
The best extinction correction for MS1512-cB58 derived from the ultraviolet continuum 
(P2000) yields E(B-V)=0.24, assuming the same extinction law; or E(B-V)=0.29 for the 
Calzetti (1994) reddening law. Both of these measurements, however, are dependent on the 
slope of the model continuum. 

We measure Ha:H / 9=3.23 ± 0.4 (note that the H/3 line is measured with SNR < 10). 
It is assumed that the underlying stellar absorption of the Balmer lines is small, and of 
equal magnitude for both lines (Olofsson 1995). Assuming Case B recombination for H I 
at T=10000 K and n e (cm- 3 )< 100, the intrinsic Ea:H.[3 ratio is 2.86 (Osterbrock 1989), 
implying E(B-V)~ 0.27 using the Calzetti law (and very little change using the LMC law). 
We measure H_7:H/3~ 0.41 ± 0.06, compared to an expected ratio of 0.45. 

The strong emission lines allow us to measure the velocity dispersion. Fitting 
Gaussian profiles to (e.g.) Ha we can measure a =FWHM/2.355. The width must then 
be corrected for the instrumental resolution. We find FWHMj^ T „ jj ne = 13.72A and 
FWHM-^ a = 20. 8A (observed). Thus, in the rest frame, a v = 81 km s -1 . Seitz et al. 
(1998) find a half-light radius ri/ 2 = lAhj^kpc. From these values, we obtain a virial mass 
M vir = 1.2 x 1O 1O M . As expected, this value is in good agreement with the masses of 
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LBGs at z ~ 3 (P98). 

3.2. Chemical Abundance 

Kobulnicky et al. (1999; hereafter KKP99) give a prescription for determining chemical 
abundance estimates from the integrated spectra of high redshift galaxies. The best 
estimates are obtained for those spectra that detect [O ill] A4363, which directly measures 
the electron temperature of the ionized gas. We do not detect this line in our spectrum, 
with a limiting flux of 0.31 x 10 _16 erg cm _2 s _1 (3cr). Reasonable (0.2dex) estimates are 
obtained in the absence of that line if the [O lll]A5007, H/3,and [O ll]A3727 lines can be 
measured, and we proceed with that limitation from the outset. 

We calculate extinction corrected ratios of [O n]/H/3= 3.52 and [O lll]A5007/H/5= 3.57. 
We then calculate the quantity R23 = (13727 + Z4959 + /5007 )/Ih(3- -R23 has a direct relationship 
to the the oxygen abundance (Pagel et al. 1979). In the case of MS1512-cB58, -R23 = 8.334. 
Using the empirical calibration of Zaritsky et al. (1994), where x = log(R2z) = 0.921, 

12 + log{0/H) = 9.265 - 0.33a; - 0.202a; 2 - 0.207a; 3 - 0.333a; 4 , (1) 

or 12 + log{0 / H) ~ 8.39. However, the relationship between i? 2 3 and oxygen abundance can 
be double-valued in this low metallicity regime (McGaugh 1991). So, following KKP99, we 
use the relationship between the [O lll]A5007/[N ll]A6584=10 116 and the oxygen abundance 
to confirm that the inferred oxygen abundance is reasonable (Edmunds & Pagel 1984). 

The solar value of 12 + log{0 / H) is 8.89, and thus cB58 has a metallicity 0.32 solar, 
which is broadly consistent with the Z/Z & ~ 0.25 found by P2000. Frayer et al. (1997) find 
a similar result, Z/Z & ~ 0.2 — 0.6, based on the spectral energy distribution (SED) and 
extinction. 

The emission line ratios also yield the nitrogen abundance in MS1512-cB58. The N/O 
ratio is typically assumed to be equal to the N + /0 + ratio. Using our estimates above, we 
can calculate N/O if we assume an electron temperature that is consistent with the oxygen 
abundance. Using KKP99's figure 5, we find T e (0 ++ ) ~ 9800 K, and then from Vila-Costas 
& Edmunds (1993; hereafter VCE), we infer T e (0+) = 11400 K from their Table 1, and 

log{N+/0+) = log F6548 + F6584 + 0.307 - 0.02 \ogt on ~ (2) 

fon ton 

where ton = T e (O+)/10 4 . Thus we find log(N + /0 + ) = -1.24. This value of the nitrogen 
abundance is well below the solar value, consistent with the metallicity inferred from the 
oxygen abundance. However, the nitrogen abundance is still higher than might be expected 
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for a galaxy as young as MS1512-cB58 has been predicted to be. The oxygen:nitrogen 
abundance ratio is most consistent with local (> lGyr old) Scd galaxies (VCE), though 
within the observed range of local, possibly younger, irregulars. Further, the relatively 
high nitrogen abundance suggests that an older population of stars is necessary to produce 
it. The ratio log(N/0) > —1.4 suggests the presence of a secondary phase of nitrogen 
production, from intermediate mass stars (Kobulnicky & Zaritsky 1999). This result 
is in conflict with spectral synthesis modeling by E96 which showed that the SED of 
MS1512-cB58 can be explained without a population of stars older than 1 Gyr. However, 
E96 only rule out an older population at the 2a confidence level, so we conclude the older 
population of stars is likely. 

3.3. Star Formation Rate 

P2000 infer an observed (without correcting for the lensing magnification factor of 30) 
SFR=2417 /if 5 2 M /yr (q = 0.1) in MS1512-cB58 from the 1500A continuum. This estimate 
included a factor of 7 correction for dust extinction (from an assumed E(B-V)=0.24). 

Kennicutt (1983) relates the SFR to Ha luminosity, assuming a Salpeter 
Initial Mass Function (IMF) with an upper mass cutoff of 100M o , by: 
SFR(M yr -1 ) = L(Ha)/1.12 x 10 41 erg s" 1 . Kennicutt assumed 1.1 magnitude 
attenuation between Ha and radio fluxes, similar to that implied by our extinction 
measurement. We neglect stellar absorption of Ha, which is likely to be only 1-5 A EW. 
We thus infer an (observed, lensing magnified) SFR=620±18 /if 5 2 M /yr (q =0.1). This 
rate is a factor of 2.18 higher than that measured by Bechtold et al. (1997), using 
narrow-band imaging. This difference is the result of an atmospheric absorption trough 
at the wavelength of redshifted Ha, that obscures approximately half the light. This 
absorption was not accounted for by Bechtold et al., who assumed z=2.7233 rather than 
z=2.729 for Ha (fortunately, their filter was wide enough to still detect the line). 

While we measure a higher SFR from Ha than Bechtold et al., we do not resolve the 
question of the discrepancy with the rate inferred from the UV continuum, a factor of 3.9 
higher. This discrepancy may be seen in the calculation of SFRs for local starbursts as 
well. Meurer et al. (1995) find UV fluxes 10 times the value predicted by Leitherer et al. 
(1995) given measured Ha; recall that similar models of Leitherer et al. (1999; the updated 
Starburst99) were also used to calculate the SFR in MS1512-cB58. Bechtold et al. speculate 
that a number of unlikely possibilities could logically account for the difference. We find no 
compelling resolution to the problem in the optical spectrum. 
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For vigorously star-forming galaxies with moderate to solar metallicity, the CO 
emission from the gas reservoir should be observable with existing ground based radio 
telescopes. However, two attempts to measure the CO(3-2) line in MS1512-cB58 have 
yielded only upper limits (Frayer et al. 1997, Nakanishi et al. 1997). The lack of detectable 
CO emission in MS1512-cB58 may be the result of either a high CO-H2 conversion factor, 
or highly efficient formation of stars from the neutral hydrogen gas. Frayer et al. suggest 
the star formation efficiency (SFE) of MS1512-cB58 is a factor of 20 larger than in local 
starbursts. However, following the arguments in Frayer et al., if we take our measurement 
of metallicity, Z/Z & = 0.32, and of SFR=62OM /yr, we would find only a factor of 1.7 
increase in the SFE. Our assumptions would not address the high L(UV)/L'(CO) ratio. We 
can, however, bring the L{Yia) / L'(CO) ratio closer to the local values (following Nakanishi 
et al.). 

4. Conclusions 

While it is difficult to draw general conclusions from observations of one object, the 
rest-frame optical spectrum of MS1512-cB58 affords us a preview of near-IR spectroscopy of 
LBGs. As expected (Seitz et al. 1998, P2000), cB58 does not appear to be a "proto-galaxy" , 
but a galaxy with significant metals produced in at least two phases. We find that the 
properties of MS1512-cB58 are consistent with the interpretation of LBGs as progenitors of 
modern-day elliptical galaxies. Previous studies have also shown that the mass, SFR, and 
extinction of LBGs are appropriate in that context (i.e. P98). We have obtained the first 
measure of metal abundance in an LBG from the promising R 2 3 method. The metallicity 
of 1/3 solar is in the range expected at this redshift from hierarchical models of elliptical 
galaxy evolution (Thomas 1999). Ellipticals produce most of their stars, and hence metals, 
in the first ~ 0.5 or 1 Gyr (Rocca-Volmerange & Fioc, 2000), and so at z ~ 3 they do not 
resemble low redshift, extremely metal poor galaxies, but rather older low-z starbursts. 
Future LBG surveys with NIRSPEC will explore this connection (Pettini et al. 2000b). 
The ease with which our spectra were obtained (1.3 hours of telescope time) demonstrates 
the great potential for observation of typical (K > 20) LBGs with NIRSPEC or similar 
instruments. 
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Table 1. Emission Lines 



Line (A) 


\ obs (/im) 




w rest \k) 


F b 


[0 n]A3726/3728 


1.1.3898 


2.72897 


37±3 


12.74±1.22 


H 7 A4340 


1.61828 


2.72875 


9±1 


1.61±0.17 


H/3 A4861 


1.81217 


2.72774 


26±4 


4.07±0.57 


[0 m] A4959 c 


1.84913 


2.72890 


26±8 


4.01±1.30 


[0 m] A5007 


1.86678 


2.72845 


97±5 


14.73±0.78 


He I A5876 


2.19100 


2.72873 


3±1 


0.35±0.09 


[0 i] A6300 


2.34949 


2.72935 


25±4 


3.06±0.46 


[Nil] A6548 


2.44204 


2.72944 


7±2 


0.86±0.21 


Ha A6563 


2.44750 


2.72935 


106±3 


12.56±0.37 


[Nil] A6583 


2.45566 


2.73031 


10±2 


1.14±0.26 



a Rest-frame equivalent width in A 

b Observed line flux in units of 10 _16 erg s _1 cm 2 

c The blue wing of the 4959A line is in a deep atmospheric absorption trough 
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Fig. 1. — N4-band spectrum of MS1512-cB58. 20 minutes of integration time. Detected 
emission lines are indicated. Dotted lines show the locations of OH sky lines with residuals 
that are visible in the 2D spectrum. The lower spectrum shows the la errors, truncated at 
14 x 10~ 17 ergs cm" 2 s _1 A _1 for clarity, and arbitrarily shifted downward on the plot. The 
increases in the errors occur at the position of the night sky lines and at positions of high 
atmospheric extinction. 

Fig. 2. — H-band spectrum of MS1512-cB58. 30 minutes of integration time. As in Figure 1, 
the lcr errors are truncated and shifted downward on the plot. The "break" in the continuum 
at 1.82 fim is a residual from the atmospheric absorption correction, not a drop in the galaxy's 
SED. 

Fig. 3. — K-band spectrum of MS1512-cB58. 30 minutes of integration time. As in Figure 
1, the la errors are shifted downward on the plot. 



